Finite-size scaling in thin Fe/Ir(100) layers 
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The critical temperature of thin Fe layers on Ir(100) is measured through Mofibauer spectroscopy 
as a function of the layer thickness. From a phenomenological finite-size scaling analysis, we find 
an effective shift exponent A = 3.15 ± 0.15, which is twice as large as the value expected from 
the conventional finite-size scaling prediction A = \jv, where v is the correlation length critical 
exponent. Taking corrections to finite-size scaling into account, we derive the effective shift exponent 
A = (1 + 2Ai)/v, where Ai describes the leading corrections to scaling. For the 3D Heisenberg 
universality class, this leads to A = 3.0 ±0.1, in agreement with the experimental data. Earlier data 
by Ambrose and Chien on the effective shift exponent in CoO films are also explained. 

PACS numbers: 05.70Jk, 64.60Fr, 68.35Rh, 75.40. Cx 
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The study of finite-size effects in critical phenomena 
has for a long time been an issue for theoretical studies, 
see Refs. |]-|3| for reviews. Only recently, due to important 
advances in epitaxial techniques necessary for the precise 
preparation of thin films, has this area become accessible 
to experiments.!-] Presently, experiments mainly study the 
finite-size shift of the critical temperature T c (n) of a thin 
film of n layers, as phcnomenologically described by the 
shift exponent A. 

Practically, two ways of defining the shift exponent 
have been used. Traditionally!!!, one measures the shift 
of T c (n) with respect to the bulk critical temperature 
T c j, = T c (oo) and sets 



ST :— (T c (oo) — T c (n)) /T c (oo) ~ n 



(1) 



in the limit when n — > oo. This defines the exponent 
A. Standard finite-size scaling theory relates this to the 
correlation length exponent v, viz. 

Alternatively, as advocated for example in Ref. |3, one 
may also define 



AT := (T c (oo) - T c (n)) /T c (n) ~ n 



-A' 



(2) 



which defines the exponent A'. Again, the limit, n — > oo is 
implied and from standard finite-size scalingut! it follows 
that A' = \jv. Empirically, using AT rather than ST 
appears advantageous, since measured values for T c (n) 
coincide with a power law for a larger range of values of 
n for AT than for ST. 

Shift exponents have been extracted successfully from 
either (|l|) or (||) for a variety of systems .ETtll For a mag- 
netic material adsorbed as a thin film on some non- 
magnetic substrate, the interaction of the magnetic mo- 
ments with the substrate may restrict the degrees of free- 
dom of the microscopic magnetic moments. Examples of 
thin films of magnetic materials or liquid heliurrusealiz- 
ing the 3r2.universality classes of the Isiiig modeH'Q, the 
XY modelQ andpfche Heisenberg modellm, as well as the 
2D Ising model,ll3 have been constructed and analysed. 



The values of the shift exponent agree with those ex- 
pected from A = A' = \ jv to within a few per cent in all 
these cases, well within the experimental error bars. In 
practice, this is an important alternative to measure true 
critical exponents, since the critical region in finite-size 
scaling studies is much larger than in thermaLrneasure- 
ments of the bulk quantities, as is well knownjjp 

However, the use of (||) has been criticized in Ref. |?], on 
the grounds that measurements qa the Neel temperatures 
in CoO/Si0 2 multilayers yieldedQ A = 1.6±0.1, but A' = 
3.4 ± 0.3, which are clearly different, ft was concluded 
in Ref. that use of eq. (||) were best avoided since the 
results for A' had no interpretation. 

In this letter, we shall re-examine this point and shall 
show how corrections to finite-size scaling may be invoked 
to provide a simple interpretation of those cases where 
a large value of A' is measured. At the same time, we 
shall report data on the critical temperature T c (n) as 
obtained from Mofibauer spectroscopy of thin layers of 
Fe on Ir(100), which can be analysed and understood 
in the same way. As a bonus, the correction-to-scaling 
exponent Ai can be measured experimentally. 

We begin describing the experimental procedure. As 
Mofibauer spectroscopy is not sufficiently sensitive to get 
some signal on so thin layers, Fe/Ir superlayers were 
prepared, where the Fe/Ir bilayer is repeated 20 times. 
In this study, the Fe thickness was varied from 2 to 
8 atomic planes and the Ir thickness is kept constant 
and equal to 1.5nm from one superlattice to another. 
The Fe/Ir(100) superlattices were prepared by molecu- 
lar beam epitaxy (MBE). All the details of the sample 
preparation are given in Ref. |l2|. The Fe/Ir system is a 
good candidate to investigate the magnetic properties of 
ultrathin films, since the Fe on Ir and the Ir on Fe growth 
modetl is an atomic layer-by-layer growth process - also 
called 2D growth. Indeed, abrupt and flat Fe/Ir and 
Ir/Fe interfaces are obtained as shown by high resolution 
microscopyEil Moreover, the occurrence of 2D growth al- 
lows us to control accurately the atomic quantity nec- 
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essary to complete one atomic layer by using electron 
diffraction. Owing to this particular behaviour, we re- 
ally deal with iron slabs constituted of exactly n atomic 
planes. A correct analysis of the variation of Fe mag- 
netic properties with the number n of deposited atomic 
planes is thus possible. In this connection, MoBbauer 
spectrometry appears well suited. This technique pro- 
vides a measurement of the so-called hyperfine interac- 
tions which occur between a resonant nucleus and its 
electronic surrounding. Among these, the magnetic hy- 
perfine interactions manifests itself by a Zeeman split- 
ting of the nuclear spin state when a steady local mag- 
netic field acts on a resonant nucleus. In case of the 
57 Fe in ferromagnetic iron, the magnetic hyperfine field 
arises mainly from the core polarization due to the 3d 
moment. The net s— spin density at the nucleus is pro- 
portional, but opposite, to the atomic on-site magnetic 
moment M. So, the room temperature MoBbauer spec- 
trum of standard bec iron consists of six line patterns; the 
measured energy positions provide the magnitude of the 
hyperfine field Bhf = 33 T. Above the magnetic order- 
ing temperature, due to fast thermal atomic spin flip in 
the paramagnetic regime, Bhf — 0: the splitting vanishes 
and the spectrum exhibits only one line (see for example 
Rcf . |l5|) . Usually, in case of normal bulk iron ferromag- 
nets, the transition takes place within a narrow tempera- 
ture range and the-thermal variation Bhf = Bhf(T) goes 
along with M(T)S3 

As only 2% of 57 Fe is present in natural Fe, the MBE 
chamber was equipped with a 92% enriched 57 Fe source, 
in order to get a sufficient amount of 57 Fe in the Fe 
layers of the superlattices. MoBbauer spectrometry was 
performed in backscattered mode by detecting conver- 
sion electron after resonant absorption of 14.4 keV 7 
rays emitted by a 57 Co source. Source drive and data 
storage were conventional, the low-temperature electron 
detector was a circular microchannel plate housed in a 
home-made cryostattll. For all investigated Fe/Ir(100) 
superlattices, the weakening and vanishing of the Zee- 
man splitting was observed to take place in a rather wide 
temperature range, compared to the collapse which oc- 
curs for standard bec iron when heating up. Thus, or- 
dering temperatures of the superlattices are here defined 
as the temperature for which the onset of a line broad- 
ening takes place in the paramagnetic spectrum when 
cooling down from room temperature. Such a determi- 
nation is illustrated in figure [l] with typical MoBbauer 
spectra recorded above and below the critical point. We 
have checked by neutron diffraction that, at least for the 
Fe thickness range where the Fe is still strained by Ir, 
there is no coupling of two Fe layers mediated the inter- 
vening Ir layer. The measured T c (n) are thus the critical 
temperatures of a single Fe layer. For each value of n, a 
single specimen was made. The quoted errors in T c (n) 
correspond to the steps in temperature used in the scan- 
ning of the MoBbauer spectra. For n = 2, . . . , 6 deposited 
iron monolayers, the ordering temperatures are given in 
tabic H For 8 monolayers, the spectrum exhibits splitting 
at room temperature. 



The dependence of the shift AT(n) on n, calculated us- 
ing the value T c ^ = 1043K for the critical point of bulk 
iron, is shown in figure 0. From at least n = 3 monolayers 
on, the data are very well described in terms of a power 
law.EJ Independently of any further theoretical interpre- 
tation, our data show that the system does lie inside the 
finite-size scaling region. Therefore, the analysis of the 
systematic variation of the critical point T c (n) with the 
number of monolayers n in terms of a phenomenological 
shift exponent is justified. We find 

A^ ff = 3.15 ±0.15 (3) 

This is about twice the value expected from A' = \jv 
in the 3Z? Heisenberg model, see table 0. Since in the 
present setting, we explore the transition from a 3D bulk 
system to a 2D film, we expect that this cross-over due 
to finite-size effects should be described in terms of the 
exponents of the 3D universality classes. 

We now come back to the problem of explaining the 
value of A' found-, rOur starting point is the theory of 
finite-size scaling.tra Given that the values of T c (n) (for 
the values of n accessible to experiment) are much lower 
than the bulk T c b, it appears sensible to consider not 
only the leading effects of finite-size scaling but to take 
finite-size corrections into account as well. For the mag- 
netization, measured on a film of a thickness of n mono- 
layers, we expecttrEl 

M n ~ n- (i l v Z (tn l ' v , un v ^ (4) 

where t ~ T c (n) — T c jj, (3 and v are the conventional 
magnetization and correlation length exponents, u stands 
for an irrelevant scaling variable which parametrizes the 
leading finite-size corrections and < is the associated 
exponent. Finally, Z = Z{z\,Z2) is a scaling function, 
which is independent of the layer thickness n. 

The pseudocritical point T c (n) (and thus an associated 
i„ as well) is determined from the vanishing of the mag- 
netization, viz. M n (t n ) = 0. This implies 

Z (t n n 1/v ,un V: ^j = (5) 

From this, we have to derive the scaling of t n as a func- 
tion of n. In particular, from figure |2| we expect that 
phenomenologically 

t n ~ t n- x (6) 

where A = A e ff is the effective shift exponent and r is a 
constant. (For the purpose of this discussion, the distinc- 
tion between A and A' as defined in eqs. (QJ|) is unnec- 
essary.) When finite-size effects are small, we can simalit 
set u = in (@) and then recover the standard resulttrEl 
A = A' = \ jv. Here, we want to discuss how large finite- 
size corrections might affect the value of A e ff. 

For n — 1 monolayers, we can find, in the same way, 
t„_i from Z(t„_i(n- l) l / v ,u(n- l) V3 ) = 0. For n large 
enough, we can expand the arguments of Z and get, to 
leading order in 1 jn 
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z 



un V3 I — y^un 



(t n n x ' u ,>, 

+ (-^n 1 /"- 1 + (t n _a - t n ) n 1 ^ Z x ~ (7) 



where = dZ / dz\^(t n n}-^ 1 ' ,un m ). To leading order in 
1/n, they can be considered as constants. The first term 
in (|?]) vanishes due to the condition (^) , used to determine 
t n . In addition, from we have t n —\ — t n ~ — Arri _A_1 . 
Inserting into (Q), we find 



A = i 
v 



2/3 



2/3 



(8) 



and, recalling that 2/3 < 0, it is already clear that the 
effective shift exponent may have values larger than the 
usually expected 1/v. The same result for A had been ob- 
tained previously from an analysis of the finite-size scal- 
ing of the correlation length.E3 

The preceding discussion can be generalized to show 
that there is a w hole scries of correction terms to finite- 
size scaling ,tL3E3 each of them with its proper value for 
the effective shift exponent A e ff = 1/v — kys, where k = 
0, 1,2,.. .. It depends on the (system-dependent and non- 
universal) value of u whether the leading correction is 
enough to describe the data or if several correction terms 
must be taken into account. For rather thick films, finite- 
size effects should be unimportant and thus k = 0, which 
reproduces the standard result A = 1/v. For thinner 
films, one might find a cross-over into a regime described 
by a non-zero value of k, where k = 1 corresponds to 
the first-order calculation given above. Higher values of 
k apply when the lower order terms vanished 

Concerning the value of w*, we now relate it to ex- 
isting theoretical predictionsEa for the three-dimensional 
0(n) model. In these models, corrections to scaling are 
described in terms of the exponent Ai, viz. M(t) ~ 
t /3 (l + at Al +...), where a is a constant. Comparing 
with the finite-size scaling form eq. (Q), then formally 
leads to [2/3 1 = A\/v. The scaling operator ip, which gen- 
erates these leading-order corrections to scaling, is well 
understood in a field theory setting of the 0(n) model. 
In particular, it is knowr£3 that rp is even under spin re- 
versal. On the other hand, the scaling operator a which 
corresponds to the magnetization is odd under spin rever- 
sal. The amplitude r of the first-order term calculated in 
(||) is proportional to the expectation value (ipcr)o, to be 
evaluated at the critical point. By symmetry, this quan- 
tity vanishes. Thus, phenomenologically, we expect the 
finite-size data of the magnetization to be determined by 
the second-order corrections, with k — 2. We therefore 
arrive at the prediction, expected to be valid in cases 
where finite-size corrections are important 



A' = A cff = ~ (1 + 2Ai) 



(9) 



In table [n[ we collect the field-theoretical predictions 
for the exponents Ai, 1/v and (1 + 2Ai)/v. The values 
quoted are the mean values of those compiled in Ref. USA 



from ^-expansion and from the resummed perturbation 
series for the three-dimensional 0{n) model. 

We now compare @ to the experimental data. First, 
consider the case of Fe/Ir(100). We expect that the fer- 
romagnetic transition of Fe is within the Heisenberg uni- 
versality class. Our result, given in eq. (|^), is in good 
agreement with the expected value A' =n3.0 ± 0.1. Sec- 
ond, we consider the CoO/Si02 system.13 CoO is known 
to be an antiferromagnet with localized moments and 
the transition is expected to be in the Ising universality 
class. The experimental result, obtained using eq. (Q), 
A c ff = 3.4 ± 0.3, is in agreement with the expectation 
A' = 3.2 ± 0.1 from table Q. Our prediction |j) is thus 
clearly confirmed for two distinct universality classes. 
To our best knowledge, this is the first time that the 
correction-to-scaling exponent Ai has been measured ex- 
perimentally. Accepting the theoretical values of v, we 
find Ai ~ 0.57(9) and 0.61(6) for the Ising and Heisen- 
berg universality classes, respectively. Our findings are 
well consistent with the theoretical predictions obtained 
from field theory, see table H 

Some more comments are in order. The consistency of 
the measured shift exponents. with A = A' = 1/v in most 
of the systems studied so faraliJ indicates that, usually, 
finite-size corrections are apparently not-.very important. 
On the other hand, in the CoO system]!] the shift expo- 
nent had been measured using both ST and AT. While 
in the second case, we have shown that finite-size cor- 
rections are needed for the proper interpretation of the 
value of A', the first case yields a value A = 1.6 ± 0.1, 
close to the expected 1/v. We stress that it depends on 
the non-universal value of the coupling u (which cannot 
be predicted from our purely phcnomcnological analy- 
sis), whether a clear power law can be observed for one 
or other of the shift eqs. (jl],^), if any, and in what regime 
of effective exponents the data will finally fall. 

In conclusion, we have presented data on the critical 
temperature of thin Fe/Ir(100) layers. The thickness- 
dependence of T c (n) has been analysed using the phe- 
nomenological theory of finite-size scaling. Finite-size 
corrections were shown to play an important role in the 
interpretation of the effective exponent value A ff, lead- 
ing to the prediction (|^), resolving an argument on the 
proper extraction of critical exponents from experimental 
data. Our data for Fe/Ir(100) and dataQ on CoO/Si0 2 
are in agreement with the theoretical predictions coming 
from the 0(n) model. This has allowed to experimentally 
confirm the field-theoretic prediction of the value of the 
correction-to-scaling exponent Ai for the 3D Ising and 
Heisenberg universality classes. 
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3 


4 


5 


6 


T c 


15 ±3 


30 ±5 


70 ± 10 


145 ± 15 


210 ± 10 



TABLE I. Transition temperatures T c (n) in Kelvin ob- 
tained from the onset of line broadening in the paramagnetic 
spectrum when cooling the system, as a function of the num- 
ber n of Fe monolayers. 



model 


n 


Aj 


l/v 


(l + 2Ai)/V 


Ising 


1 


0.50(2) 


1.586(4) 


3.17(6) 


XY 


2 


0.53(2) 


1.492(7) 


3.07(11) 


Heisenberg 


3 


0.555(25) 


1.413(9) 


2.98(13) 



TABLE II. Field theory predictions^! from the three-di- 
mensional 0(n) vector model of some critical exponents. The 
numbers in brackets estimate the error on the last given digits. 
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FIG. 1. Example of the Mofibauer line broadening when 
the temperature is decreased, for a 6ML 57 Fe superlattice. 
Down to T = T c , the FWHM F is constant. When T c is 
reached, F begins to increase. 
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FIG. 2. Shift AT from eq. (|2j) of the critical temperature 
of the Fe/Ir(100) system as a function of the number n of Fe 
monolayers, for n = 2, . . . , 6. The curve is the power-law fit 
AT(n) = 1069 n -315 . 
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